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ABSTRACT 

We compare the isotropic equivalent 15 — 2000 keV 7-ray energy, i? 7 , emitted by a sample of 91 swift 
Gamma-Ray Bursts (GRBs) with known redshifts, with the isotropic equivalent fireball energy, Er,, 
as estimated within the fireball model framework from X-ray afterglow observations of these bursts. 

The uncertainty in E 1 , which spans the range of ~ 10 51 erg to ~ 10 53,5 erg, is « 25% on average, due 
mainly to the extrapolation from the BAT detector band to the 15 — 2000 keV band. The uncertainty 
in Afb is approximately a factor of 2, due mainly to the X-ray measurements’ scatter. We find E 1 and 
25fb to be tightly correlated. The average(std) of 7/ lhr = log 10 (15 7 /(3£ e Elf b lhr )) are —0.34(0.60), and 
the upper limit on the intrinsic spread of ry 7 is approximately 0.5 (e e is the fraction of shocked plasma 
energy carried by electrons and E^f* is inferred from the X-ray flux at x hours). If the uncertainties 
in the determinations of and Eft, are twice larger than we estimated, then the data imply no 
intrinsic variance in ?y 7 . We also find that Er> inferred from X-ray observations at 3 and 11 hours are 
similar, with an average(std) of log 10 (U^ lr /l?f b lhr ) of 0.04(0.28). The small variance of % implies that 
burst-to-burst variations in s e and in the efficiency of fireball energy conversion to 7-rays are small, 
and suggests that both are of order unity. The small variance of ?y 7 and the similarity of E^ r and 
£’f b lhr further imply that s e does not vary significantly with shock Lorentz factor, and that for most 
bursts the modification of fireball energy during the afterglow phase, by processes such as radiative 
losses or extended duration energy injection, are not significant. Finally, our results imply that if 
fireballs are indeed jets, then the jet opening angle satisfies 9 > 0.1 for most cases. Extending our 
analysis to late times we find a significant reduction in £*,, < E^/E^ >= 1.4, consistent with jet 
breaks on a 1 d tim e scale in a significant fra ction of the bursts. These results are consistent with the 
main results of lFreedman fc Waxmanl (120011) . which were based on a much smaller sample of GRBs. 
Subject headings: cosmology: observations- 7-ray: sources- 7-ray: bursts 


1. INTRODUCTION 

GRBs are the most powerful explosions in the Uni¬ 
verse, and include the highest redshift objects observed. 
The widely accepted phenomenological interpretation of 
these co smological sources is the so called ’’Fireball (FB) 
model” (iPaczviiskilllDSfil [Goodmanlll986lh In’’optically 
thin” versions of this model, the energy carried by the 
hadrons in a relativistic expanding wind (fireball) is dis¬ 
sipated through internal shocks between different parts 
of plasma. These shocks reconvert a substantial part of 
the kinetic energy to internal energy, which is then radi¬ 
ated as 7-rays by synchrotron and inv erse-Compto n ra¬ 
diation of sh ock-accelerated e lectro ns ([Rees fc Meszarosl 
Il992t iNaravan. Paczvnski fe Piranll 1 9921) . Alternatively, 
photospheric models for the prompt emission have b een 
proposed (e.g. iPe’er et all l2012t iBeloborodovI I2013D as 
being significant contributors to shaping the prompt 
emission. Regardless of the mechanism responsible for 
the emission of the prompt gamma-rays, the fireball 
is expected to drive a shockwave into its surrounding 
medium, which decelerates as it encompasses an increas¬ 
ing amount of mass and is believed to be responsible 
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for the afterglow emission on time scales of minutes to 


years (iPaczvns 

ri & Rhoads! 19931: iMeszaros & Rees 1997t 

IWaxman 1994 

Sari et al. 1998h. Lone term afterglow ob- 


servations led to the concl usion (Rhoads 1999: lSari et al.l 
Il999t iRacusin et al.l [20091) that the fireball is not spheri¬ 
cal bur rather jet-like. 


The model for the afterglow phase is completely de¬ 
fined by the fireball energy, Er,, the jet opening angle 
0, the fraction of shocked plasma energy carried by elec¬ 
trons and magnetic field, e e and £b , the spectral index 
of the power-law energy distribution of shock accelerated 
electrons, p, and the density of the plasma into which 
the shock expands, n. While afterglow observations are 
in general in good agreement with the model, it remains 
difficult to tightly constrain the parameters of the model 
based on observations. This is due to the fact that a 
complete calibration of the parameters requires a deter¬ 
mination of the three breaks in the spectrum (cooling, 
peak and self absorption frequencies), which in turn re¬ 
quires a multi wavelength coverage of the afterglow over 
hundreds of days. In particular, the fireball energy, and 
hence the efficiency with which this energy is converted to 
7-rays, can be determined in only a few cases, for which 
adequate spectra l coverage is available (e.g. GRB970503, 
iFrail et al.ll2000fl. 

iFreedman fc Waxmanl (120011 . hereafter FW01) have 
shown that it is possible to estimate the fireball energy 
(note that throughout this work, ’fireball energy’, as well 
as GRB energy, refer to the isotropic equivalent energies) 
carried by electrons (e e Efb) using a measurement at a 
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single time of the flux at a frequency for which the emis¬ 
sion is dominated by fast cooling electrons (i.e. electrons 
for which the cooling time is shorter than the dynamical 
expansion time), and that this estimation is weakly de¬ 
pendent on poorly constrained model parameters. The 
physical basis underlying this estimate is that the lumi¬ 
nosity at such a frequency corresponds to the rate at 
which energy is deposited in shock accelerated electrons, 
~ e e Efh/t. Since such observations have been carried 
out for many GRBs, this allows one to infer the effi¬ 
ciency of 7 -ray production for a large sample of GRBs. 
Using this method, it was found in FW01 that E 1 and 
e e Efb, inferred from X-ray observations at ~ 11 hours, 
are strongly correlated, implying that the burst-to-burst 
variations in e e and in the efficiency of fireball energy 
conversion to 7 -rays are small, and suggesting that both 
are of order unity. It was also pointed out that if GRB 
fireballs are indeed jet-like, then the jet opening angle 
should satisfy 9 > 0.1 for most bursts (for smaller open¬ 
ing angles, sideways expansion of the jet would occur 
at earlier times, thus significantly affecting the inferred 
isotropic equivalent energy, at t < 11 hours). While these 
results are significant, they were based on a small sam¬ 
ple of GRBs and their afterglows: a total of 13 GRBs of 
which only 7 had measured redslrifts. 

The applicability of the results and conclusions of 
FW01 to the majority of GRBs is challenged by more re¬ 
cent studies, which suggest that the fireball energy is sig¬ 
nificantly modified on a time scale of hours following the 
burst by e xtended energ y injection or significant radia¬ 
tive losses. iNousek et alJ (l2006ll suggested, based on an 
analysis of 9 lightcurves from a sample of 27 GRBs, that 
the fireball energy is increased by a factor > 4 due to en¬ 
ergy injection on a timescale mostly up to 3 ho urs (with 
two exception of up to 5 and 11 hours); IZhang et all 
(200(f) argue that a good fraction of GRBs (though they 
do not quantify how many) have shallow decay phases 
on timescales of up to a few hours, and interpret them 
as corresponding to continual energy injection, which in- 
creases the fire ball ener gy by a factor of up to 10 , and 
iPanaitescu fo Vestrandl (120121 ) find significant energy in¬ 
jection to be ubiquitous in afterglows through analyz¬ 
ing the consistency of breaks with the fire ball model. 
Analysing two G RB afterglows iBerger et al.l (120041) con¬ 
cluded, based on lYost et al.l (120031) .that 50% and 90% of 
the fireball energy of the two bursts was radiated away 
during the early (3 hou rs and 6 d ays) afterglow. It was 
furthermore argued (jSari fc Esinl 12001 11 that the X-ray 
flux based estimate of the fireball energy Er, should be 
strongly affected by inverse-Compton (IC) energy loss of 
the electrons, which may lead to a significant fraction of 
the electron energy being radiated at (unobserved) fre¬ 
quencies well above the X-ray band (where the emission 
is dominated by synchrotron radiation), and which is ex¬ 
pected to strongly vary between bursts. 

Such significant modifications of the fireball energy (or 
of the estimated fireball energy, e.g. due to IC suppres¬ 
sion of the synchrotron emission), which vary from burst 
to burst and are not expected to be correlated with the 
prompt 7 -ray emission, would introduce a large scatter 
to the Ery/Efb ratio, which would be inconsistent with 
the results of FW01. In this paper we expand the anal¬ 
ysis of FW01 to a large sample of Swift GRBs, in order 
to examine whether the tight correlation of C 7 and An,, 


and the implications of such a tight correlation, hold for 
the majority of the GRBs. The large sample and the im¬ 
proved afterglow data available enable us to quantify the 
uncertainties more reliably and thus also to draw quan¬ 
titative conclusions. 

The GRB sample used in our analysis is described in 
§[2] The methods used for estimating E 1 and Er>, and 
the uncertainties in these estimates, are described in § [3] 
and in § [4] respectively. The correlation between i ? 7 and 
Ufb is analyzed in § [5] The implications of the relatively 
tight correlation between E 1 and Er, are discussed in §E1 
and our conclusions are summarized in § [3 

2 . THE GRB SAMPLE 


Our sample consists of 91 long GRBs with known red- 
shifts, which were detected by Swift in the period be¬ 
tween May 2006 and August 2009 (from GRB060502A 
to GRB090812). The total number of long GRBs that 
were detected in this period, and for which redshifts are 
known, is 112 (7 short, Tgo < 1.8 s, GRBs with known 
redshift that were detected during this period are not in¬ 
cluded in our sample). 21 of the 112 GRBs were excluded 
from the analysis due to the following reasons: for 1 GRB 
the required BAT data are not available; 14 GRBs do 
not have satisfactory XRT data at the 3-11 hours time 
interval (For 2 of these bursts no XRT data are avail¬ 
able, for 8 bursts no data are available at t > 4 hr, and 
for 4 there are only 3 data points in the time range of 
3hr < t < 11 hr. We note that the lack of data for 
these bursts is unlikely to be due to low X-ray fluxes, 
and is likely due to low sampling cadence, since the ex¬ 
isting data points show relatively high fluxes.); 6 bursts 
have late time (close to 3 hours) flares or clearly rising 
lightcurves around the relevant times. 

The exclusion of the latter group of 6 GRBs from our 
analysis implies that our conclusions may not apply to a 
small minority, ~ 5%, of the long GRB population. 


3. ESTIMATING E FB 

The fireball energy carried by electrons can be esti¬ 
mated from a measurement of the flux density f v at a 
frequency v, which is above the cooling frequency (corre¬ 
sponding to the frequency of radiation emitted by elec¬ 
trons with cooling time equal to the dynamical time), 
using equations 4 and 5 of FW01: 

Ee-Ufb = (c 2 c' 3 )- 1/2 c 1 - 1 -^^/ 1 ,(^,t)y E , (1) 

1 + Z 


where 




£ = 


P~ 2 
p + 2 
( 2 ) 


and Ci,( 72,(73 are numerical constants taken from 
FW01, Ci = 1.4 x 10 ~ 21 cm 3 / 2 , C 2 = 6.1 x 

10~ 5 s 3 / 2 g - 1 / 2 cm -1 (for p = 2.2) and C 3 = 6.9 x 
10 39 s _ 3 /, 2 g _ 1 ,/ 2 cm _2 . The inferred energy is independent 
of the density of the plasma into which the fireball ex¬ 
pands, and very weakly dependent on the value of £b 
for p « 2. For our nominal energy estimates we use 
£ e — £b = 1/3. The dependence on £ e is shown explic¬ 
itly in all our results. 

For the spectral index, we adopt in our analysis a value 
of p = 2.2. We show below that such a universal value 
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is supported by both the afterglow spectra (§ m and 
the time dependence of afterglow X-ray flux (§ E2). As 
discussed in some detail in § 16.31 the tight correlation 
between Er, and E 1 further supports a universal value 
of p. Using a value of p = 2.4 instead of p = 2.2 increases 
the estimate of Er> by 60% (FW01). 

Eqs. ([Tim yield an accurate estimate of £ e -Ufb under the 
assumptions that the emitting plasma is well described 
by the self-similar spherical fireball model and that elec¬ 
tron acceleration is well characterized by time indepen¬ 
dent e e and p. When these assumptions hold, the X-ray 
flux should drop with time as fx oc i( 2-3 U/ 4 , and Er, 
inferred using eq. © should be independent oft. A mea¬ 
sure of the accuracy of the determination of e e Efb, and 
a test of the validity of the above assumptions, is there¬ 
fore obtained by comparing the values of £ e Efc, inferred 
from the X-ray flux at significantly different times. When 
the assumptions described above are not valid, eqs. m 
provide a less accurate estimate of £ e Efb. 

We show in § 13.21 that the values of E^ in¬ 
ferred using eq. © at 3 hr and If hr are simi¬ 
lar (implying dlogfx/dlogt ~ — 1 ): we find that < 
log 10 (i?^ lr /£ ; ^ b lhr ) >= 0.04 and that only « 20% of 
the bursts have E^ r and £j! b lhr values differing by more 
than a factor of 2. These results support the validity of 
the model assumptions described above. They may ap¬ 
pear to be contradictory to the claims that ’’flat” X-ray 
lightcurves (with d log fx/d log t significantly larger than 
- 1 ), and he nce signific ant energy injection, are common 
Ce.g. lNousek et al’ll2006il . This issue is discussed in § 13.31 

Finally , we note that the uncertainty in the energy es¬ 
timate should include a systematic uncertainty due to the 
uncertainty in the values of the constants {Ci, C 2 , C 3 }, 
which vary by factors of a few between different models 
of the afterglow synchrotron e missio n (see for example 
the discussion in lGranot. Konigl fe Piranll2006f) . We do 
not include in our analysis a treatment of this system¬ 
atic uncertainty, since a modification of {Ci, C 2 , C 3 } will 
modify the values of Efb and of E 1 /E{ b by a fixed factor 
for all bursts, but will not affect the fractional scatter of 
these values. It will thus not affect the main conclusions 
of this work. 

3.1. Constraints on p from afterglow spectra 

Following FW01, we present below evidence for the 
validity of the approximation of p ss 2 based on the spec¬ 
tral index of the afterglow emission, determined from the 
X-ray and optical fluxes, /3 ox = -Hf x /f 0 )/ln(v x /v 0 ). 
At times when the cooling frequency is below the X-ray 
band, we expect (jp — l )/2 < ft ox < p/2 (with the lower 
limit obtained when the cooling frequency is in the X-ray 
band, and the upper limit when it is in the optical band). 

In fig. |T] we show f3 ox at 11 hours for the first 38 
GRBs in the samp le (up to GRB071031), as given by 
iZheng et all (1200911 . For 75% of the bursts, /3 ox is in the 
range of 0 . 6 - 1 . 1 , which is the range expected for p = 2 . 2 . 
Of these bursts for which the measured spectral index 
does not lie in this range, one has f3 ox > 1.1 while the 
others have /3 ox < 0.6. Such values are in principle not 
consistent with the shape of the spectrum in the fire¬ 
ball model for p = 2.2 at times when the cooling fre¬ 
quency lies below the measured X-ray band. However, 
as already noted by FW01, dust extinction in the host 
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Fig. 1.— A histogram of the effective spectra l index /3 ox fo r a sub - 
set of 38 bursts within our sample, given bv IZheng et al.l (|20090 . 

galaxy, which suppress es the optical flux, can reduce the 
observed value of j3 ox . IZheng et al.l (|2009h find that low 
values of f) ox typically correspond to high values of hy¬ 
drogen column density, which serves as a proxy for the 
strength of local dust extinction. Moreover, they state 
(referring to previous work) that for several of these low 
j3 ox bursts, detailed studies of the optical afterglow al¬ 
lowed one to measure the optical extinction values, thus 
raising the value of j3 ox . Finally, we have compared the 
group of bursts with low f3 ox in the sample to the whole 
sample and found both groups to have the same average 
.Efb and E^/Er,, thus indicating that the the low fi ox 
bursts do not represent a separate population. 

3.2. fx and Er, at 11 and 3 hrs 

For each GRB we have estimated the 10 keV flux den¬ 
sity at 11 and 3 hours after the GRB trigger by fit¬ 
ting a power-law temporal behavior to the Swift X-ray 
lightcurve between 3 x 10 3 s and 6 x 10 4 s. Due to the 
rapid decay and to the flares often characterizing the X- 
ray light curve at early time, we do not use t < 3 x 10 3 s 
data (except for 4 GRBs, which clearly show a regular 
behavior of the lightcurves at earlier times and for which 
the earlier time data are necessary in order to obtain a 
reliable interpolation at 3 hours). Due to the possible 
appearance of a jet-break suppression of the flux at late 
time, t > 1 d, we do not use t > 6 x 10 4 s data (except 
for 12 GRBs, for which a late break in the lightcurve is 
clearly absent and for which later data are necessary to 
obtain a reliable interpolation at 11 hours). 

The uncertainties in the derived X-ray fluxes at 3 and 
11 hours are dominated by the scatter in the reported X- 
ray flux measurements, as illustrated in figure [ 5 ] We use 
this scatter as indicative of the flux uncertainty, rather 
than the reported flux uncertainties, since the latter do 
not appear to be consistent with the scatter of the data 
points. The fluxes at 3 and 11 hours are estimated us¬ 
ing a least square fitting of a first degree polynomial in 
the logt — log fx plane, and their uncertainties are de¬ 
termined such that 90% of the data points lie within the 
uncertainty. We note that for 20% of the bursts, for 
which the last measured X-ray data point is at a time 
earlier than 5 hours after the trigger, the extrapolation 
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t [sec] 

Fig. 2. — XRT flux measurements for GRB090618. The straight 
line is the power-law time decay fit for the range of interest, based 
on which the fluxes at 3 and at 11 hours are inferred (red circles). 
The red lines represent the estimated flux uncertainties (see text). 
The empty circles represent data points at t > 6 x 10 4 s, which are 
not taken into account in determining the fluxes at 3 and 11 hours. 



Fig. 3.— The relation between E fb inferred from the X-ray fluxes 
at 3 and 11 hrs. The dashed (solid) blue lines correspond to the 
energy estimates at 3 and 11 hrs being within a factor of 2 (10) of 
each other. 


of the power law fit to 11 hours makes a significant con¬ 
tribution to the uncertainty. 

The average(std) of d log fx/d log t we obtain from our 
fits to the X-ray light curves between 3 x 10 3 s and 6 x 
10 4 s are -1.25 (0.52), consistent with the value of -1.1 
expected for p = 2.2. The deviation of the decay rate 
from this value is larger than the estimated error in the 
determination of the decay rate for < 20 % of the bursts. 
This result is illustrated in figure [3] showing the tight 
correlation between the fireball energies estimated at 3 
and 11 hours. The average(std) of log 10 (.E^ r /-Ef b lhr ) are 
0.04(0.28). Given our estimated uncertainties in inferring 
Efbi the data are consistent with no intrinsic variance in 
the ratio E^ r /E^ hT . Out of the 91 bursts in the sample, 
11(14) have a fireball energy estimated at 3 hours which 
is less than half (more than twice) that estimated at 11 
hours. We show in § [5] that the average of 7? 7 is similar 
for bursts with small (less than a factor of 2 ) and large 
(more than a factor of 2 ) differences between E and 
Efo hr , while the variance of t ? 7 is somewhat larger for the 
latter group. 


3.3. Flat light curves 

The results reported in the preceding paragraph may 
appear to be contradictory to the claims that ’’flat”, 
i.e. d\ogfx/d\ogt significantly larger than -1, X-ray 
lightcurves (a nd h ence significan t energy injection) are 
common (e.g. INousek etal.ll2006f). Th is is, h owever, not 
(necessarily) the case. INousek et al.1 (2006) concluded 
that the canonical lightcurve contains a shallow decay 
phase based on an analysis of 27 bursts out of which only 
9 showed flat decay phases, and only 5 (< 20%) showed 
such behavior extending to t > 6 x 10 3 s. Our analysis, 
based on a much larger and complete sample, shows that 
the fraction of afterglows showing a flat phase leading to 
a factor 2 or larger increase in the inferred E fb between 
3 and 11 hr is < 10%. 

At times earlier than 3 x 10 3 s, the lightucrves are 
in general less regular, as was discussed in § 13.21 with 
many of the lightcurves showing both steeper and flatter 
time dependence compared to later times. In particular, 
roughly 25% of the bursts show lightcurves which are sig¬ 
nificantly flatter at early times compared to late times. 
However, as we explain in detail in § 16.61 the contribu¬ 
tion of possible early energy injection to the total fireball 
energy cannot in general be significant. 


4. ESTIMATING E 1 

The 7 -ray spectrum of_GRBs is typci ally described us¬ 
ing a Band function ([Band et al.lll993l) . 


dN 

~dE 


E(2 + c) 

T.50 KeV^ 6 *' Peak ^E < -Ebreak ( 3 ) 

i®>£break 


where Ab rea k = ^ ( 9 +a) eak ■ Due to the limited 

energy range of the BAT detector, 15 to 150 keV 
(! Bart helm v et al.1 [2005 1. the typical values inferred for 
the Band function parameters based on BAT data dif¬ 
fer from the typical values inferred based on data from 
BATSE, which is sen sitive in the range of 30 to 2000 keV 
llKaneko et al.I [2006 b While the average low energy in¬ 
dex a in the BAT sample used here, —1.11, is simi¬ 
lar to the average value of the BATSE GRBs, —1.08 
(jKaneko et al.l l2006i i, the peak energy A pea k and the 
high energy index /3 inferred from the BAT and from 
the BATSE data are quite different. In particular, the 
average peak energy obtain ed for BATSE spectra is 
260 keV (|Kaneko et al.ll2006l h outside the BAT energy 
band. Thus, in order to determine E 1 in the energy 
range of 15 to 2000 keV for the BAT GRBs, we use a 
Band function with a and flux normalization as given 
in the Swift data archive (determined by the BAT spec¬ 
trum), but use the average values of the B ATSE catalog , 
Ep e ak = 260 keV and /3 = —2.4 (jKaneko et al.l 12006 ). 
This procedure introduces of course an uncertainty in 
the estimated value of A 7 , due to the uncertainty in the 
extrapolation of the spectrum up to 2000 keV. We esti¬ 
mate the uncertainties using extrapolations to high en¬ 
ergy with the extreme values obtained in the BATSE 
catalog for A pea k and /3: 150 keV and 500 keV for Apeak 
and —2.6 and —2.0 for /3. 

The method for estimating E 1 and its uncertainty de¬ 
pends on the values of the Band function parameters 
given in the Swift data archive, as explained in detail 
below. We note that in all cases, the low end of the 
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model spectrum we adopt is identical to that found by 
BAT (and its flux integrated over the 15-150 keV range 
matches the BAT reported integrated flux). 

• For /3 in the range of [-4,-2], which holds for « 15% 
of the bursts in our sample, the flux per logarithmic 
photon energy interval, E 2 dN/dE, is decreasing 
with E. We therefore estimate the [15,2000] keV 
fluence by integrating the BAT Band function up 
to 2000 keV, and the uncertainty range by extrap¬ 
olating the BAT spectrum beyond 150 keV using a 
/? = —2.6 and a j3 = —2.0 power-laws. 

• For /? > —2, which holds for Ri 55% of the bursts, 
E 2 dN/dE is increasing with E. We therefore es¬ 
timate the [15,2000] keV fluence by extrapolating 
the BAT Band function to 260 keV, followed by a 
j3 = —2.4 power-law at higher energy. The lowest 
value of energy uncertainty range is obtained by 
extrapolating the BAT spectrum beyond 150 keV 
using a /3 = —2.6 power-law, and the highest value 
by extrapolating the BAT Band function to 500 
keV followed by a /3 = —2.0 power-law at higher 
energy. 

• For /? < —4, which holds for « 30% of the 
bursts, inbreak inferred by the SWIFT analysis is 
well above the BAT band. We therefore estimate 
the [15,2000] keV fluence by extrapolating the BAT 
Band function to 260 keV, followed by a /3 = —2.4 
power-law at higher energy. The lowest value of the 
energy uncertainty range is obtained by extrapolat¬ 
ing the BAT Band function to 500 keV followed by 
a /? = — 2.6 power-law, and the highest value by 
extrapolating the BAT spectrum beyond 150 keV 
using a /3 = —2.0 power-law. 

The ratios of the [15,2000] keV fluences, inferred as de¬ 
scribed above, and the BAT [15 150] keV fluences are in 
the range of 1.02 to 4.26, with an average of 2.25. The 
uncertainties due to the extrapolation are ~ 25% on av¬ 
erage, reaching a factor of 2 for some bursts. In addition 
to the uncertainty introduced by the extrapolation to 
high energies, we include in our analysis the uncertainty 
in the BAT fluence as reported in the SWIFT archive. 
The average value of the latter uncertainty is 7%, and 
it amounts to no more than 26% for any burst in the 
sample. 

5. THE CORRELATION BETWEEN E 1 AND E FB 

Using the methods described in the preceding sub¬ 
sections we find the average(std) of log 10 (£’ 7 /lerg) and 
of log 10 (3e e £ , f b lhr /lerg) to be 52.3(0.7) and 52.5(0.6) re¬ 
spectively, very close to the values inferred in FW01. As 
can be seen in figure[4] log(U 7 ) and log(-E’fb) are linearly 
correlated, with a high correlation coefficient of rs 0 . 6 . 
For the average(std) of ? 7 7 lhr = log 10 (E 1 /(3£ e E^ hl )) we 
find —0.34(0.60), consistent with the results of FW01 
0.01(0.5) (and also with those of iD’Avanzo et all 120121 
who give their nominal results with a normalization of 
£ e = 0.1). 

Given the estimated uncertainties in the derived values 
of E-f and of e e E^ hr , the hypothesis that their ratio is 
universal (i.e. the hypothesis of no intrinsic variance in 


TABLE 1 

Average values and Standard deviations 


time 

set^ 

sample size 

3 hr 
all 

91 

3 hr 
tight 
66 

11 hr 
all 

91 

11 hr 
tight 
66 

< logloU^A-Efb) > 

52.6 

52.6 

52.5 

52.6 

var(log 10 ((3e e )Efb)) 

0.6 

0.6 

0.6 

0.5 

<V~t> 

-0.36 

-0.43 

-0.34 

-0.41 

var(?j 7 ) 

0.55 

0.51 

0.60 

0.51 

X 2 pdf^ 

3.3 

2.9 

3.6 

2.6 

intrinsic variance^ 

0.41 

0.36 

0.50 

0.36 


a ”all” = the entire burst sample analyzed in this work, ” tight” = 
a subset for which the fireball energy estimates at 3 hours and at 
11 hours are within a factor of 2 of each other. 
b X 2 , as defined in the text. 
c required to obtain x 2 = 1. 

? 7 7 lhr ) is inconsistent with the data (resulting in % 2 per 
degree of freedom of ~ 4). In order to estimate the intrin¬ 
sic variance of ? 7 7 lhr , we assume that it follows a gaussian 
distribution with a variance s. s is determined by equat¬ 
ing x 2 t° 1 in the equation y 2 = (N — 1 )~ 1 S y *~^~ 2 T?7 , 
where Xi,yi are log(U 7 ), log(iffb), and eq represents the 
uncertainty in the energy estimates (of = a 2 . + of.) 
(minimizing y 2 while setting s = 0 gives the nominal 
value of ?) 7 , which is then used in computing s). We find 
the intrinsic variance to be approximately 0.5. We note, 
though, that this value is sensitive to the uncertainties 
in Ej and Er,. For example, if these uncertainties are 
twice larger than we estimated, then the data implies no 
intrinsic variance in ?y 7 . 



Fig. 4.— The relation between E-y and The GRBs marked 

with filled circles are those for which the fireball energy estimates 
at 3 hours and 11 hours lie within a factor 2 of each other. We 
show for comparison the results obtained by FW01 (red diamonds). 
Solid lines correspond to F 7 /(3c, ; )/f) lir =0.1,1.10 and the dashed 
line to -E 7 /l?(3£ e )-E 7b lhr = 0.01. 

Repeating the above analysis using X-ray obser¬ 
vations at 3 hours, we find the average(std) of 
log 10 (3e e F;^ hr /lerg) and ?yf hr = log 10 (U 7 /3 e e E^) to be 
52.6(0.6) and —0.36(0.55) respectively, similar to the val¬ 
ues estimated at 11 hours. 

Table [T] gives the values of the average and variance of 
the fireball energy and of ?y 7 , obtained with and without 
bursts for which the fireball energy estimates at 3 and 11 
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hours differ by more than a factor of 2 . As can be seen 
from the table, the results obtained with and without 
these bursts are very similar. We note that, as shown in 
Fig- HI the variance of ? 7 7 is larger for bursts characterized 
by larger differences between E^ r and E^ hr (while its 
average value is nearly independent of this difference, 
Fig. [7|. 



Fig. 5.— Same as fig. [4] for i?fb inferred from the X-ray flux at 
3 hrs (instead of at 11 hrs). 


While rj-y lies within the range of ~ 0.1 to ~ 1 for 
most of the GRBs in our the sample, there are 5 GRBs 
(GRB080607, GRB080319B, GRB061110B, GRB061007, 
GRB060510B) with exceptionally large values of ?y 7 , 
77 7 ( 3£ e ) -1 > 8 . Out of these, 4 bursts have exception¬ 
ally steep lightcurves between 3 and 11 hours relative to 
what is expected from the fireball model. One possible 
explanation for such behavior could be that these GRBs 
suffer from significant radiative losses at these times, im¬ 
plying that their E-y/Efo is actually lower than inferred 
by using the estimates of E^ T and Uf b lhr (Estimating the 
radiative losses for particular bursts must be done using 
extensive lightcurve fitting to constrain the fireball model 
parameters, which is beyond the scope of this work). It 
should however be pointed out that, as shown in Fig. 0 
for these bursts the ratio of our inferred [15,2000] keV flu- 
ence to the [15,150] keV measured BAT fluence is 4, the 
largest in our sample. This implies that the E 7 and ? 7 7 
estimates of these bursts are more sensitive to systematic 
errors in our extrapolation method. 

6 . DISCUSSION 

As mentioned in the introduction, significant modifi¬ 
cations of the fireball energy, or of the estimated fire¬ 
ball energy, which vary from burst to burst and are not 
expected to be correlated with the prompt 7 -ray emis¬ 
sion, would introduce a large scatter to the E 1 /E^ ratio. 
Thus, the relatively tight correlation we find (§ 0 be¬ 
tween U 7 and Eft, sets relatively tight constraints on the 
combined effect of processes that lead to such modifica¬ 
tions. Quantitatively, the intrinsic variance in ?y 7 , which 
is inferred from the scatter in the E^ - E^ correlation, 
sets an upper limit of a factor of 2-3 on the scatter of the 
Ej/Efb ratio due to the combined effects of these pro¬ 
cesses. In what follows we discuss the main implications 
of the tight correlations we find between E^ and Ef b and 



Fig. 6.- The standard deviation of 77-y of that part of the bursts 
within the sample for which the fireball energy estimated at 3 and 
11 hours is within a factor of 10 x . 



Fig. 7. — The average value of rj<y of that part of the bursts within 
the sample for which the fireball energy estimated at 3 hrs and at 
11 hrs are within a factor of 1CP of each other. 

between E^ hr and E^ T . 



Fig. 8. — versus the fraction of the total estimated GRB en¬ 
ergy which is contained in the BAT range. Empty circles represent 
bursts for which E^ hr and differ by more than a factor of 

2. The solid line shows the average of 7y 7 , the dashed lines are 
separated from the average by the variance of 


6.1. Jet opening angle and ”off-axis” detections 
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The fact that the fireball energy estimates at 3 and 11 
hours are similar indicates that in general jet breaks do 
not occur well before 11 hours, since that would cause 
the inferred E^ hl to be on average significantly lower 
than E^ r . Since the relation between the jet opening 
angle 9 , and the jet br eak time tj is estimated as (e.g. 
iLivio fc WaxmaiJfeOQOf l 

6^ 0.12 (E isOt53 /n 0 )- 1/8 t 3 /J ay , (4) 

we conclude that typical opening angles are greater than 
0.1, as inferred in a similar manner by FW01. 

A similar argument also implies that the population of 
observed GRBs can not be dominated by jets observed 
initially (during y-ray emission) ”off-axis”, i.e. with an 
angle 9 between the line of sight and the jet axis, which 
is larger than 9j. In such cases, the evolution of the fire¬ 
ball and of the emission of the afterglow radiation are 
not described by the simple self-similar model leading to 
eqs. CLED, which are used to infer E{b, the luminosity 
of the prompt y-ra y emissi on would vary strong ly de¬ 
pending on 9 (e.g. iLundman Pe’er fc Rvdel 1201311 . and 
the tight correlation between E 1 and Ef b would be a co¬ 
incidence. 

In order to test for the possible presence of jet breaks at 
t > 11 hr, and to demonstrate the sensitivity of our anal¬ 
ysis method to such breaks, we have derived E ^ for the 
73 bursts within the sample for which X-ray light curves 
extending to t > 2 d (with at least 2 measurement points 
at t > 1 d) are available, using data points at t > 1 d only. 
For these bursts we find < log 10 (£^ lr /.E^ i ) >= 0.16, 
significantly larger than the uncertainty in estimating 

< log 10 Ea, >, which is « 0.07. The hypothesis that 
the average energies at 3 hr and 2 d are equal is re¬ 
jected at a ~ 85% confidence level based on a standard 
two sample t-test (Recall that in §[3]we have found that 

< log 10 (£’^ u '/£’^ i ) >= 0.04, consistent with no change 
in < log 10 Efb > between 3 hr and 11 hr). A possible ex¬ 
planation of the « 40% reduction in < E fb > betweenn 
3 hr and and 2 d is a jet break at t < 2 d. A precise 
characterization of the nature of jet breaks is beyond the 
scope of this work, but the « 40% drop could correspond 
(in two extreme cases) to a universal jet break time of 
« 30 hr or to ss 1/2 the bursts having a jet break at 
11 hr. 

6.2. Variations of the "f-ray emission on small angular 

scales 

Due to the relativistic expansion of the fireball, the 
flux we observe is obtained from a conical section of the 
fireball (around the line of sight) with an opening angle of 
1/T, where T is the (time dependent) expansion Lorentz 
factor. The GRB y-rays are expected to be emitted at 
the highly relativistic, r = T 7 ~ 10 2 5 , phase of fireball 
expansion, while the X-rays observed at ~ 10 hrs are 
expected to be emitted after deceleration to T == Tx ~ 
10 1 . Thus, the isotropic equivalent fireball energy, Efb, is 
determined by an average of the jet properties over angles 
~ l/Tx ~ 0.1, much larger than those from which y-ray 
emission is observed, ~ l/r 7 < 0.01. If the emission 
of y-rays were to vary significantly over angular scales 
0.01 ~ l/r 7 < A 9 < l/Yx ~ 0.1, a large scatter would 
be introduced to the E^/Er, ratio. The absence of such 
large scatter implies that the jet’s y-ray emission does 


not vary significantly over 0.01 < A 9 < 0.1 angular scales 
within the jet opening angle 9j. 

6.3. Variability in the microscopic parameters of the 
afterglow model: £ e , p 

Significant burst to burst variations in e e and/or p 
would lead to a significant scatter in the E- f jE^ ratio, 
while significant dependence of £ e and/or p on the shock 
Lorentz factor would lead to a significant deviation from 
unity of the ratio E^ hl /E^ r . Our results therefore in¬ 
dicate that e e and p are uniform between bursts, which 
is reasonable given the fact that they are determined by 
the micro-physics, and that their values do not strongly 
depend on the shock Lorentz factor. 

We note that some earlier stu dies have found signifi¬ 
cant scatter in e e . For example, iYost et ah ([2003h find 
values in the range of 0.12 — 0.34 by fitting multiband 
observ ations to 4 different afterglows, and iBerger et all 
(120041) find a factor of 4 difference between the two af¬ 
terglows for which they present their analysis. While the 
intrinsic variance we find does not exclude such scatter 
in £ e , it would imply that almost all the intrinsic vari¬ 
ance we find in the f ? 7 / Eft, ratio is due to the variance in 
£ e . This, in turn, would imply that the variance due to 
burst-to-burst variations in y-ray production efficiency 
is negligible. Since £ e is a parameter determined by mi¬ 
croscopic processes, while the y-ray production efficiency 
depends on the macroscopic properties of the flow (e.g. 
Lorentz factor variability within the fireball wind, see 
SlOi. this appears to be unlikely. 

6.4. Variability in the GRB j-ray production efficiency 

Strong burst-to-burst variations in the efficiency of the 
conversion of the fireball energy to y-rays would lead to 
a large scatter in the E 1 /Ep a ratio. The absence of such 
scatter implies that the efficiency is not highly variable 
between bursts. This is a strong constraint on the fireball 
model, since significant variations are natural to expect 
in most of its variants. For example, in models where 
the emission of y-rays follows from the conversion of fire¬ 
ball kinetic energy to internal energy by internal colli¬ 
sions within the wind, the efficiency strongly depends on 
the amplitude and structure of Lorentz factor variations 
within the expanding fireball wind. Uniform efficiency 
suggests an efficient process (with order unity efficiency) 
of converting fireball energy to radiation. 

6.5. Radiative losses 

The tight correlation between E^f hl and E^f *, and 
the fact that the ratio E^ hT /Ef^ T is consistent with 
unity, implies that no significant radiative cooling oc¬ 
curs between 3 and 11 hrs. In addition, the small in¬ 
trinsic variance in E^/Efb implies that if radiative losses 
are significant at earlier times, they are nearly uniform 
among all GRBs. These observations can be used to 
impose constrains on parameters of the fireball after¬ 
glow mode l. W e dem onstrate this by using the model 
of IBerger et alJ (120041) for radiative losses, in which the 
fireball energy drops with time as £*, oc t ~ 17£ / 12 with 
£ = £ e f (1 + 1.05e e ), as long as all electrons are in the 
fast cooling regime. 

We first note that radiative losses cannot account 
for most of the scatter between fireball energy esti¬ 
mates at 3 and 11 hours, since the average and std 
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of log 10 (£ l fb l 3hr/-E l fb,iihr) are 0.04 and 0.28 respectively, 
while radiative losses can only make this ratio larger than 
unity. Thus, we assume that only the average value is af¬ 
fected by radiative losses, i.e. on average no more than 
20% of the fireball energy is radiated away between 3 and 
11 hours. This implies either that £ e is small, e e < 0.15, 
in which case radiative losses are small altogether, or 
that the stage at which all electrons are in the fast cool¬ 
in g regime ends at t < 1 hr hour. Using the estimate 
of iSari et al.l (|1998f ) for the duration of this fast cooling 
stage, ffast, together with our estimate of £ e Efb, we find 
n = 2 x 10“ 2 (£ B /0.1) _2 (£e/0.3) _1 (ff ast /l hr)cm -3 . 

To conclude, the small difference (< 20%) between 
the estimated fireball energies at 3 and 11 hours im¬ 
plies either e e < 0.15, or a circum burst density n < 
2 x 10 -2 (£ B /0.1) -2 cm -3 . 

6.6. Energy injection 

The tight correlation between -Ey b lhr and E^j r , and the 
fact that the ratio E^ hr /E^ r is consistent with unity, 
implies that no significant injection of energy takes place 
during this time period. Can significant energy injection 
at yet earlier time be a common feature of GRBs? 

Since some the afterglow lightcurves contain periods 
of shallower decay than expected by the basic fireball 
model, it has been suggested that additional energy 
might be injected to the fireball at times later than the 
initial burst. Based on their analysis of 9 lightcurves 
(out of a samp le of 27) containing a phase of shallow 
decay. iNousek et all (120061) suggest that energy injection 
might increase the fireball energy by a factor > 4 on a 
timescale mostly up to ss 3 hours. In fact, their analy¬ 
sis (see their table 3) suggests that for these burst, the 
increase in fireball energy can even be much more signif¬ 
icant, on the order of a factor of 10 on average (which, if 
taking into account radiative losses which occur simul¬ 
taneously, could actually correspond to an even more 
significant energy injection). The variance in the fac¬ 
tor b y w hi ch the e nergy is increased is ~ 2 — 6. Simi¬ 
larly, [Zhang^Fal] (120061) present lightcurves containing 
a shallow decay phase, and interpret them as indications 
for late time energy injection, but they do not quantify 
how representative of the ” average GRB” these examples 
are. In a m ore recent w ork focusing on energy injection, 
iPanaitescu fc Vestrandl 12 0121 1 find that among a sample 
of ss 100 swift afterglows featuring a break in the X-ray 
lightcurve, « 30% of the breaks can be explained as be¬ 
ing jet breaks in an adiabatic model, while « 60% can 
be explained as jet breaks in a model including extended 
energy injection, thus suggesting that energy injection on 
a time scale of up to ~ 1 day is a ubiquitous feature of 
GRB afterglows. 

Our results, on the other hands, favor the possibility 
that late time energy injection is not a significant fea¬ 
ture in GRBs. First, the relatively low intrinsic variance 
which we find in the ratio of the prompt energy and of the 
fireball energy, E^/Er, implies that any significant en¬ 
ergy injection must be correlated with the prompt emis¬ 
sion, i.e. it should increase the fireball energy by a con¬ 
stant factor for all bursts. While we cannot rule out such 
a possibility, it would mean that the already challenging 
high efficiency of conversion of fireball energy to 7-rays 
should be even higher, which would then be difficult to 



10 3 10 4 10 5 

t [sec] 

Fig. 9.— XRT data for GRB090516. The straight line is the 
power law fit, from which the flux values at 3 and 11 hours are 
inferred (red circles). The uncertainties in the inferred values (red 
lines) are based on the scatter of the X-ray data points. Open cir¬ 
cles denote data points at times which are outside the time window 
which we use for determining the X-ray fluxes at 3 hr and 11 hr 
(3 x 10 3 s to 6 x 10 4 s for all bursts in the sample). 

explain. Moreover. IPanaitescu fc Vestrandl (j2012f ) find a 
large scatter, of « 2 orders of magnitude, in the factor 
by which energy injection increases the fireball energy, 
which is inconsistent with the small variance we find in 

Ej/Efb. 

Seco nd, th e significant ene rgy injection inferred by 
IPanaitescu fo Vestrandl ( 20121 : during the time period of 
100 s to 100 kilo-s implies an average increase in energy 
by a factor of 1.5 — 3.5 between 3 hours and 11 hours, 
inconsistent with our finding that E^ hr and E^ T are the 
same on average. Third, their analysis is based on de¬ 
tecting breaks in the lightcurves, which are interpreted as 
jet breaks which occur at t < 11 hr for 60% of the cases. 
This is inconsistent with our finding that jet breaks do 
not typically occur at such early time. 

Finally, looking specifically at t he bursts co mmon to 
our sa mple and to the sample of IPanaitescu fe Vestrandl 
(I2012T) . and which feature a > 3 hours break, which ac¬ 
cording to their analysis can only be explained with the 
presence of late time energy injection, we find the esti¬ 
mated fireball energies at 3 and 11 hours to be similar. 
As an exam ple, t he X RT data for GRB090516, for which 
IPanaitescu fc Vestrandl (120121) find a break at 2 x 10 4 s, 
is shown in figure [9] In our analysis, the energy esti¬ 
mates at 3 and 11 hours are similar for this GRB, and 
the break, which is not analysed in the current work, 
appears to occur later. 

6.7. Inverse-Compton losses 

If the electrons responsible for the emission of X-rays 
via synchrotron emission lose a significant fraction of 
their energy by inverse-Compton (IC) emission at much 
higher frequencies, this would lead to a significant under 
estim ate o f the fireball energy Er. As was already noted 
by I Berger et al.l (f2003l ), since the effect of IC losses de¬ 
pends strongly on n (and on £ B ), we would expect it to 
vary strongly between bursts, leading to a large scatter 
in Ej/Er,. The small scatter in E^/Er suggests that IC 
losses are not significant at 3 and 11 hrs. We do not carry 
a detailed analysis of the implications of this conclusion 
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to fireball model parameters, but note that IC emission 
may be suppressed by the Kl ein-Nis hin a effect for reason¬ 
able model parameters fe.g. iNakar. Ando fe Sari!l2009l) . 


7. CONCLUSIONS 

We have analyzed a sample of 91 swift Gamma-Ray 
Bursts (GRBs) with known redshifts, and found E 7 
and Ufb to be tightly correlated and E^ r and Uf b lhr 
to be tightly correlated (§ 0. The average(std) of 
7y 7 lhr = log 10 (U 7 /(3£ e Lyb lhl )) are —0.34(0.60), and the 
upper limit on the intrinsic spread of r]^ is approximately 
0.5. If the uncertainties in the determinations of £/ and 
Ufb are twice larger than we estimated, then the data 
imply no intrinsic variance in r\ 1 . The average(std) of 
log 10 (£'^ lr /£ , ^ ) lhr ) are 0.04(0.28). Given our estimated 
uncertainties in inferring Er,, the data are consistent 
with no intrinsic variance in the ratio E^f 1 /E^ 1 . 

The implications of these results were discussed in de¬ 
tail in § 0 The small variance of ?? 7 implies that burst- 
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